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Abstract 

We compute logarithmic corrections to the 0(a 2 ) lattice artifacts for a class of lattice 
actions for the nondinear O(n) sigma- model in two dimensions. The generic leading 
artifacts are of the form a 2 [ln(a 2 )]^ n ^ n ~ 2 ^ . We also compute the next-todeading 
corrections and show that for the case n = 3 the resulting expressions describe well 
the lattice artifacts in the step scaling function, which are in a large range of the 
cutoff apparently of the form 0(a). An analogous computation should, if technically 
possible, accompany any precision measurements in lattice QCD. 



1. Most of our knowledge concerning renormalization of quantum field theories stems 
from perturbation theory. Although there are no rigorous proofs in general, many 
of the results are structural and hence considered to carry over to non-perturbative 
formulations. Indeed there is supporting evidence from various studies, e.g. of 
soluble models in 2 dimensions and of 1/n expansions of some theories. 

The same situation holds concerning cutoff artifacts in lattice regularized the- 
ories. It is generally accepted that these artifacts are summarized in Symanzik's 
effective action. In this framework generic lattice artifacts are, in particular for 
asymptotically free (or trivial) theories, expected to be integer powers in the lat- 
tice spacing 0(a p ) ,p = 1, 2, . . . up to possible multiplicative logarithmic corrections. 
This is an extremely important ansatz in the extrapolation of lattice data to the con- 
tinuum limit, especially for present computations of lattice QCD where the lattice 
spacings are typically around O.lfm. 

In this letter we will examine in more detail Symanzik's theory for the 2- 
dimensional non- linear O(n) a- model. In such a simple bosonic model lattice ar- 
tifacts are expected to be of the form 0(a 2 ). It was thus rather surprising that 
precision measurements [T], some years ago, of certain observables in the 0(3) sigma 
model exhibited apparently linearly dependent 0(a) artifacts for a rather large range 
of computable lattice spacings. The importance of finding the solution of this puzzle 
was emphasized by Hasenfratz in his lattice plenary talk in 2001 [2]. 

To define the measured quantity referred to above, one considers the model 
confined to a finite (1-dimensional) box of extension L (with periodic boundary 
conditions). The LWW coupling [3] is defined as 



where m(L) is the mass gap of the theory in finite volume. Next one measures u\, 
defined similarly with doubled box size. In the continuum limit u% is a function of 
no, called the step scaling function u\ = cr(2, uq). For the lattice regularized theory 
there are lattice artifacts and 



The advantage of this measurement for the purpose of studying lattice artifacts is 
that there is no need to know the box size L or the mass gap m(L) in physical units. 

The results of the MC measurements are shown in Fig. [TJ One can see that 
the lattice artifacts (cutoff effects) are very nearly linear as function of the lattice 
spacing a both for the case of the standard lattice action (ST) and for a modified 
action (MOD). Although the effects are in this case relatively very small, they seem 
not of the theoretically expected form. Note however the encouraging feature that 
computations with different lattice actions are consistent with the same continuum 
limit, supporting the crucial concept of universality. 



uq = L m(L) 



(1) 



ui = 2L m(2L) = S(2, u , a/L) . 



(2) 



1 



1 



.31 



W 1.27 



o 



1.29 




1.25 







0.05 



0.1 



a/L 



Figure 1: Monte Carlo measurements of the step scaling function at uq — 1.0595. The data for 
the larger artifacts correspond to a modified action. The fit contains a and a In a terms. 

The 2-dimensional 0(3) model is integrable and the finite volume mass gap 
(and hence the step scaling function) is exactly calculable using thermodynamic 
Bethe Ansatz techniques [3]. However, it turned out that the knowledge of the 
exact continuum limit did not help much in clarifying the problem of artifacts. This 
is shown in Fig. [2J where the exact continuum values are already subtracted. The 
"linear" fit assuming the functional form 5(a) = c±a + C2alna + c^a 2 , still gives 
a much better representation (x 2 /dof = 0.9) than a "quadratic" fit of the form 
5(a) = c\a 2 + C2a 2 \na + c^a 4 which has an unacceptable x 2 /dof = 7.7. Note that 
although the L/a > 10 data are used in the fits, the "linear" fit describes well also 
the three coarsest points. 

In the early 80's Symanzik was working on the nature of lattice artifacts, in 
particular with respect to his improvement program [5 6 7 8]. In this letter the 
general theory is not discussed; we only consider Symanzik's theory applied to the 
2-dimensional O(n) a- model. Nevertheless the spirit of the general theory can al- 
ready be understood by studying this example. We shall see, based on Symanzik's 
theory, why in the 2-dimensional cr-model quadratic artifacts are expected and how 
the particular logarithmic corrections predicted by this theory solves the puzzle of 
apparent linear artifacts. Similar conclusions have been reached previously in stud- 
ies of the O(n) model in the first orders of the 1/n expansion |9|10] . We can here 
only give a brief description; full details will be presented in a separate paper [TT] . 
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Figure 2: MC measurements of the step scaling function at no = 1.0595. The curves are "linear" 
(dashed line) and "quadratic" (solid line) fits described in the text. 

2. We write the lattice Lagrangian including the source terms symbolically as 

Aatt = i ( 9 » S ■ ^) latt -J-S, S 2 = l, (3) 

with some lattice regularization of the kinetic term. This is used in the generating 
functional for bare lattice S'-field correlation functions, which, after Fourier trans- 
formation, become functions of the momenta, the bare lattice coupling Ao and the 
lattice spacing a. Performing, for fixed momenta and coupling, a small a expansion 
we can write them as 

g i x att (X ,a) = g x ^(X ,a)+a 2 g x ^(X ,a) + O(a 4 ) , (4) 

where the upper index X symbolizes any r-point correlation function (in x-space 
or in Fourier space), and both the scaling functions Q x (°) and the leading cutoff 
corrections Q x ^> are still weakly (logarithmically) depending on a. 

The separation of the full lattice correlation function into a scaling piece and 
cutoff corrections is unambiguous and straightforward in perturbation theory (PT). 
In fact, in PT at I loop order both terms are finite (order £) polynomials in In a. 
One of our main assumptions here is that the expansion makes sense also beyond 
PT. Usual renormalization theory deals with the scaling part Q x \°) . Symanzik's 
important contribution was to show that the next term Q x ^ can be generated 
using an effective Lagrangian. 
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£ = —(d fi S-d tl S)--I-S, (6) 



3. Symanzik's local effective Lagrangian C e s can be specified in the continuum in 
D = 2 — e dimensions in the framework of dimensional regularization: 

7 

- £ cff = -£ + a 2 !■(<?, £)[/,, (5) 
i=i 

where C is the continuum Lagrangian with source terms 

—2 (d^S ■ d^S) - -g. 

and the U{ form a basis of local operators of dimension 4 invariant under the lattice 
symmetries, discussed below. There are no such operators of dimension 3 and hence 
no 0(a) terms. 

First we recall the continuum part ([6]). Here go is the bare coupling and the 
square of the bare 0(n) field S a (x) is normalized to unity which is parameterized, as 
usual, by S l = goir 1 , i = 1, . . . , n — 1; S n = a = yjl — <7q7t 2 . The source dependent 
action A. = J d D xC(x) is used in the generating functional 

Z[I] = J (Vn) e- 4 , (7) 

which can be used to obtain bare correlation functions of the field S a (x): 

e—c, ...,«,)- ^'M^ ■ ■ -jj^l w . («) 

where the functional derivative is taken at I§ (x) = m^S 0,11 ; i.e. a mass term (external 
magnetic field) is introduced to avoid infrared singularities. As is well known, for 
0(n) invariants the limit mo — > can be taken at the end of the calculation. 

In their seminal paper [i!2] Brezin, Zinn-Justin and Le Guillou prove the renor- 
malizability of the 0(n) model using functional methods. They showed that the 
generating functional 2 _1 [Io]^[^] is finite as a function of the renormalized quanti- 
ties j a (x),g, fj,, mji if we write 

I a (x) = Z 1 (g,e)Z- 1 / 2 (g,e)g 2 f(x)+lS, 

g 2 =f, e Z 1 (g,£)g 2 , (9) 
m 2 = Z 1 (g,s)Z- 1 / 2 (g,e)m 2 R , 

where in the minimal subtraction scheme the renormalization constants contain only 
pole terms in e. 

Functional derivation with respect to the source j a (x) gives renormalized corre- 
lation functions, i.e. correlation functions of the renormalized fields S% = Z~ 1 / 2 S a : 

g? m (g,l*,e) = Z- r / 2 (g,e)g x (go,e). (10) 
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We assume that X is 0(n) invariant and that the mR — > limit has been taken. 
Finiteness means that the limit 

G(r)(9,v) = ^9f R) {g,^e) (11) 

exists and defines the renormalized correlation function in two dimensions. 

The renormalization group (RG) equations express the fact that the bare cor- 
relation functions are independent of the renormalization scale /i. In terms of the 
renormalized correlation functions this is expressed as 

{v + ~j(g)}gf R) (g,tJ,) = 0, (12) 

where the RG differential operator is 

p = 4 +,3<9) l' (13> 

and the RG beta and gamma functions are defined as 

P(g) = f - - £ d iz l{9 ,e) = -/W - Pig 5 -P2 9 7 + ... (14) 
1 2 + g — df 2 ^ 

( \ I n( \ e $\ dlnZ (a, £ ) 2,4, / 1K x 

i{g) = \P{g) - Y j — ^ — = 7o5f +715 +■■■• ( 15 ) 

Next, the dimension four operators U{ in (0) are linear combinations of the 5 
Lorentz scalar operators considered by Brezin et al |12j : 



o 1 = I (d^s ■ d^s) 2 , o 2 = l (d^s ■ d u s) (d^s ■ d u s) 



o 3 = ~ns-ns, o 4 = \ad^s ■ d^s , o 5 = \a 2 

where 



(16) 



a = , (17) 

a 

and a further two operators which are only invariant under discrete lattice rotations: 

D D 

(i=i fi=i 
where t and k are the traceless parts of the symmetric tensors t, k: 

t/xupa = S • d^dvdpdaS , 

1 (19) 
k^pa = - {(d^S ■ d u S) (d p S ■ d a S) + 2 perms} . 
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Although the source dependent operators O4 and O5 look O(n) non- invariant, Brezin 
et al show that they must be included in the operator renormalization scheme for 
consistency. The operators renormalize multiplicatively according to 

U i{R) =JCij(g,e)Uj, (20) 

where the matrix of renormalization constants is block diagonal, consisting of a 2 x 2 
block for i,j = 6, 7 (oc A, B), and a 5 x 5 block for i,j = 1,2,3,4,5. The operator 
renormalization matrix is of the form 

2 4 4 

ICij (9, e) = Sij ~ 9 — kij + y £ "if + ^2 ( k ^sj + 2^o%) + ■ • • • (21) 

Our one loop result kij for the 5x5 sub-block agrees with that which is obtained 
from the computation of Brezin et al |12j in the basis Q16|) . 

Renormalized correlation functions with one operator insertion are given by 

g? {R) (g lt i) = limZ- r / 2 (g, £ )JC ij (g,e)Gf(go,e) . (22) 
They satisfy the RG equation 

{^ + ^(9)}Gf(R)(9^) + ^j(9)Gf iR) (9,^ = 0, (23) 
where the anomalous dimension matrix is defined by 

"M = IC ts (9,e) (/3(g) - f ) ^ ^ ^ = -hrf+v^ g±+v\f gH. . . . (24) 

We have chosen the basis 11% such that the one-loop anomalous dimension matrix 
kij is diagonal of the form 

hj = 2(3 Ai5ij , (25) 

with eigenvalues corresponding to 

f n 1 — n 1 ] , „ 

l n — 2 n — 2 n — 2 J 

Now defining 

7 

Cj(g,s) = Y J ^(9,e)ICT j 1 (g,e) (27) 

8=1 

we can rewrite the second term in ([5]) as 

7 7 

Y J YiU i = Y,CiU i{R) . (28) 

8=1 8=1 
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The limit Cj(<?, 0) must therefore exist: 

oo 

c t (g)=c t (g,0)=Y,cf ) 9 2e . (29) 
i=o 

The coefficients q depend on the particular lattice action. In our investigations 
we only considered actions quadratic in the spins: 

A =\ EE^Ws-yWy), (30) 

x,y a 

with j3 = 1/Ag. Here K is short range, satisfying ^2 x K(x) = 0, and K(z) = 
K{Rz) , where R is a lattice rotation or reflection. For the standard action K(z) = 

E M [ 2 ^,o - 8 Xt fi - 5 z _p] . 

The tree level effective action for this class of lattice actions involves only the 
operators A and O3: 

- ct s - -555 < 8 - s • 8 - s » + f IS- 4 + >} + (° 4 ) ■ < 31 > 

where eo = = 1 for the standard action, and off-shell improved actions are such 
that eo = e4 = 0. From this we obtain directly the leading coefficients c[°^ in (|29j) . 
4. With these preparations the precise relation between the lattice correlation func- 
tions and those obtained by using the effective action can now be specified as 

g x ^(X ,a) = y r (g)Gf R) (g,^j , 

i=i v a 

where we have identified (for simplicity) the scale parameter /i of dimensional regu- 
larization with the inverse of the lattice spacing a. The finite wave function renor- 
malization constant y(g) comes from the relation j a (x) = y{g)J a {x) between the 
lattice source J and the (renormalized) dimensional regularization source j. 
The result (f32j) can also be written as 

Ctt(Ao,a) = e x(0) (Ao,a) {l + a 2 5 x (X , a)} + O (a 4 ) , (33) 

where 

7 

5 x (X ,a) = Y J C l (g)Sf(g,a) (34) 

i=l 
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xXf \ ^i(-R) ( S ' a) ,ok\ 

h (9, a) = x 1 n • (35) 
Gfn) (9, a) 

It is easy to see that the functions 5 X satisfy the RG equation 

{~ a lL + (3{9) l^} S i C (9,a) = -"iMSf(g,a). (36) 

To derive an explicit expression for S x we must solve this partial differential equa- 
tion. For this purpose we introduce the matrix Uij(g), which solves the ordinary 
differential equation 

Ui j (g) = -p is (g)U sj (g), (37) 

where 

'~ iKgJ ~~ f-^ j ' 

If we find the solution of (|37p we can write the general solution of (|36p as 

8f(g,a)=U lj (g)Df(A), (39) 

and the lattice artifacts are of the form 

7 7 
5 x (X ,a) =J2*i(9)D?(A), Vi(g)=J2 c MU si (g). (40) 

i=l s=X 

The functions D x depend only on A, the RG invariant combination of g and o. These 
functions are non-perturbative and depend on the quantity (X) we are considering. 
On the other hand, the coefficients Vi(g) are perturbative and they remain the same 
for all physical quantities (but depend on the lattice action we started with). 
We take the following ansatz: 



u ij (g) = \Si j + T / ^9 2e - 2 \ 9 - 2A \ (41) 

i=2 



Here the coefficients kf^ may still weakly (logarithmically) depend on the coupling. 
This can arise if the difference between two eigenvalues Aj — Aj is a non-zero integer, 
which is possible in our case, i.e. for n = 3. We will however ignore this subtlety, 
because we have verified that for the quantities we need here it plays no role. 
We can now write the lattice artifacts as 

7 

5 x (X ,a) = Y,<9)g- 2A > D?(X) , (42) 
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where 

oo oo 

Vl (g) = Ci {g) + £ c s ( 5 ) J]) *2 } ^ 2 = £ ^ ^ • ( 43 ) 
s 1=2 e=o 

The spectrum of one-loop eigenvalues given by (|26p plays a crucial role in our con- 
siderations. The leading term corresponds to 

Tl 

Ai = = n X = 1 + 2 X , (44) 

n — 2 

(where x := l/( n ~~ 2)) and the subleading one to 

1 



X- (45) 



n - 2 

We thus have the leading expansion 

5 x (\ ,a)= Vl D? (g- 2 ) 1+2x + v 5 DX (g-y + ..., 

(46) 



^{( 5 - 2 ) 1+2x ^ 0) +^ (<r 2 ) 2x } + o(( 5 



-2\X\ 



It turns out that for the n=3 case Vg = Cg = 0. This means that in this case the 
corrections start one power later and we have the leading expansion 

^(A , a) = D? {,i 0) <T 6 + v? g^ + 5 " 2 } + 0(1) . (47) 

The first expansion coefficients are 



„(o) _ jo) _ e o 

Vl ~ Cl ~4(n-l)' 



Concretely we have 



(48) 



- i „(2) m 

sl - 2(A 1 -l-A s ) Ps1 ' { yj 

which is different from zero for s = 1, 2 only. 

This and the connection between the lattice coupling Ao and g is all we need to 
write down the final result: 

**<*•■■> ^ 4 ( „-i)Vr ^< A > { (*) 1+ * + ^ } + (*) (») 

for n > 4, and 

<^(Ao, a) = ^_ ^ (A ) j> + r (2)^2 + r ( 3 )^j + 0(1) ^ (51) 
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for n = 3, where we have introduced the inverse coupling j3 = 27t/Aq. We expect 
this form of artifacts to be generically present for all observables. We note that our 
final result, eq. (|50p . is completely consistent with the large n results of ref. [9]. 
We have computed the coefficient which is composed of three contributions: 

r (2) =r (2) +r (2) +r (2)^ (52) 

For n = 3 it would be nice to know the 3-loop coefficient ; its computation would 
however be a major undertaking 0. 

For the computation of the first term in ()52[) 



rf = , (53) 

we need the 1-loop coefficients of the effective action. We obtained these by calcu- 
lating the 2- and 4-point functions in both regularizations. More precisely we need 
the scaling part and the 0(a 2 ) piece of the lattice correlation functions, and in the 
continuum we need the original correlation functions as well as the ones where those 
dimension four operators that appear in the tree level effective action are inserted. 
This is a long computational, the details of which will be published in ref. [llj . 
The second term in (|52p 

^ = ^(1-V0, (54) 

involves the 1-loop relation between the lattice coupling Ao and the renormalized 
coupling of dimensional regularization, g: 

g 2 = Xl + ^\t + ... (55) 

This has been known for a wide class of actions for a long time. e.g. for the case of 
ST at n = 3 one gets rf) = -3.98028. 

Finally the third term can be obtained from the 5x5 two-loop anomalous 
dimension matrix of the dimensionally regularized scalar operators 



n-2 11 n 



(56) 



I It is built from (among other things) the two-loop coefficient appearing in the effective 
action and the three- loop anomalous dimension matrix elements vf- . 

2 In the case of operator insertions there are many Feynman diagrams to be calculated. In the 
4-point function case for example there are 11 non-trivial diagrams, which cannot be reduced to 
simpler ones like renormalization of tree diagrams or insertion of 2-point function subgraphs with 
operator insertion. Each of these is a complicated function of the four momenta and since they are 
topologically distinct, it is difficult to automate the calculation. 
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This is again a lengthy computation which leads to the simple result 

(57) 



.(2) _ o 9 



111 2(n-2)' 

5. For the standard action we finally obtain 

/n\ 5.6416 n , . 

r ( 2 ) = -0.7625 -, (58) 

n — 2 4 

giving a very large negative coefficient, —7.1541, for n = 3. Using the result (|5ip . 
we can write in this case the leading terms of the asymptotic series describing the 
lattice artifacts in terms of a and the inverse coupling: 

const, a 2 [P 3 - 1.1386/3 2 + 0(/3)] . (59) 

Due to the big negative value of the subleading term, the coupling depen- 
dent function in the square bracket is a very rapidly growing function, which can 
compensate the decreasing of one power of the lattice spacing, resulting in a be- 
havior approximately linear in a in this limited range of coupling between (3 ~ 1.6 
and (3 ~ 2.0 as observed in pQ. To demonstrate this, we multiplied the measured 
lattice artifacts by L 2 to remove the a 2 factor and fitted a polynomial cubic in 
(3 to the resulting function. We first made a two-parameter fit, with the result 
3.60 [/? 3 — 1.52/3 2 ] . We also tried a 2-parameter fit, where the relative coefficient of 
the subleading term was fixed at the value we obtained from our analysis; we found 
3.06 [(3 3 - l.U/3 2 - 0.60/?] . This fit is shown in Fig. El For both fits the two small- 
est lattices were omitted. Both fits are good, indicating that the MC measurements 
are consistent with the functional form given by our result (|59p . 

For the modified action MOD, where in addition to the nearest neighbor cou- 
pling there is also coupling between spins in the diagonal direction (with equal 
strength), the results are similar. For this lattice eo = 5/3 and 

^ 7.2625 3n , , 

A 2 > = -3.0153 + — , (60) 

ra- 2 20 v ; 

giving —9.8278 for n = 3. The artifacts are given by a formula similar to (|59p . 
where the overall constant is bigger by a factor 5/3. For the relative (MOD/ST) 
artifacts the overall non-perturbative constant cancels and we find the following 
parameter-free 1-loop prediction (for n = 3): 



MOD/ST= - 



1 + ^ + 0(1//^ 



(61) 



which varies between 1.85 and 1.89 in our limited range of coupling. This is in very 
good agreement with the ratio of artifacts found in [lj, as can be seen in Fig. [H 
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Figure 3: Modified 2-parameter fit: 3.06[,3 3 - 1.14/3 2 - 0.60/3]. The points are from lattice sizes 
L/a = 5, 6, 10, 12, 16, 24, 32, 64; the fit excludes the first two. 



6. We conclude that Symanzik's theory describes well lattice artifacts in the 0(3) 
sigma-model. By looking into the details of this problem we had the opportunity to 
further gain confidence in Symanzik's theory of artifacts and improvement. Similar 
computations should, in our opinion, accompany precision lattice studies of QCD in 
order to control and better estimate systematic errors arising from lattice artifacts 
for extrapolations to the continuum limit. 
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